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1. Introduction 

Many important aspects of the low energy dynamics of QCD are related to the chiral symmetry 
and its spontaneous breaking in the QCD vacuum. Lattice QCD is a promising tool to actually 
simulate the low energy regime of QCD and to reproduce the spontaneous symmetry breaking and 
related phenomena, provided that the chiral symmetry is preserved on the lattice. This was the 
prime motivation for the JLQCD and TWQCD collaborations to embark on the QCD simulations 
with exact chiral symmetry, i.e. dynamical overlap fermion simulations, despite the enormous 
computational costs they require. 

Theoretically, with the dynamical overlap fermions, one may investigate the relation between 
chiral symmetry breaking and low-lying modes of the Dirac operator, as suggested by the Bank- 
Casher relation and the chiral random matrix theory. Also, the quantities related to the topology of 
the SU(3) gauge field, such as the topological susceptibility, may be studied using fermionic proves 
using the singlet axial-Ward-Takahashi identity. 

Phenomenologically, the chiral symmetry plays an important role in the chiral extrapolation 
of lattice data for many physical quantities. Since the use of the continuum chiral perturbation 
theory is justified, the chiral extrapolation is largely simplified when the chiral symmetry is exactly 
preserved on the lattice. In this talk, I will discuss a few examples that demonstrate the chiral ex- 
trapolation performed using the next-to-next-to-leading order chiral perturbation theory formulae. 
We have other new applications for which the exact chiral symmetry plays a crucial role. They in- 
clude a calculation of the VV —AA vacuum polarization functions and an extraction of the strange 
quark content from the (partially quenched) nucleon masses. 

This project was started in 2006 when the present supercomputer system, that comprises Hi- 
tachi SRI 1000 and IBM Blue Gene/L, was installed at KEK. A report of the project at an earlier 
stage was presented by Hideo Matsufuru at Lattice 2007 which mainly described the simula- 
tion itself. The present talk discusses about the physics applications from the project in some detail. 
For individual contributions, see their own write-ups [Q, ||, ^, |5|, |^ ^ |8|]. 

2. Simulation strategy and status 

The overlap-Dirac operator D(0) = p/a[l+X/ VX^X] (with the Wilson kernel X = aDw — p) 
1^, 1^] provides a theoretically ideal fermion formulation on the lattice, as it has an exact chiral 
symmetry at finite lattice spacings while satisfying the index theorem. Implementation of the 
overlap operator is however non-trivial, since the definition contains a discontinuity at the zero of 
the denominator X^X (or the zero of the hermitian Wilson-Dirac operator H-w = Y5X). One usually 
approximates the sign function sgn(//vi') using polynomial or rational functions, but the near-zero 
eigenmodes must be specially treated. At finite lattice spacings, it is known that there are finite 
density of near-zero modes of H-w jr2|], that means that the computational cost to identify the near- 
zero modes grows as quickly as C?(V^) and prevents one from performing large scale simulations 
using the overlap fermion. 

The near-zero modes are associated with a dislocation of the background gauge field and thus 
are unphysical (an explicit example is found in [pj[]). We suppress them by introducing extra heavy 
Wilson fermions (and associated ghosts) [|l^, [l^]. They produce a factor det[//^/(//^ + jx^)] to the 
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Boltzmann factor (/i is a mass of the associated ghosts) and naturally suppress the near-zero modes. 
The numerical cost for the overlap fermion simulation is then dramatically reduced especially when 
the lattice volume is large. 

Since the topology change (or the change of the number of zero-modes of the overlap-Dirac 
operator) may occur only when a small eigenvalue of changes its sign, the suppression factor 
strictly prohibits the tunneling among different topological sectors, as far as the Monte Carlo algo- 
rithm is based on a continuous change of the gauge field configuration. This is in fact a property 
of the continuum QCD, i.e. topology change cannot occur through continuous deformation of the 
gauge field. 

Dynamical overlap fermions have been attempted by several authors |T^, 17, [l^], in which the 
treatment of the topology tunneling is one of the main issues. Among them, our simulations have 
reached the level to produce broad physics results for the first time. On a lattice of size ~ (2 fm)^ we 
simulate dynamical fermions with several quark masses between ms/6 and with the physical 
strange quark mass. In addition to the runs for two-flavor QCD, which is already published [jl^], 
we have completed a series of runs for 2-1-1 -flavor QCD. The main reason that made this possible 
is the topology-fixing, as well as the computational power provided by the BlueGene/L. 

The two-flavor runs were done at j8 = 2.30 using the Iwasaki gauge action on a 16^ x 32 lattice. 
The lattice spacing determined through the Sommer scale ro = 0.49 fm is 0.118(2) fm. For each 
of six sea quark masses covering the region 171^/6 ~ m^, we accumulated 10,000 HMC trajectories 
in the trivial topological sector Q = 0. (For a check of the fixed topology effect, we also produced 
configurations of Q = —2 and —4 at one sea quark mass. For details, see below.) 

The 2-1-1 -flavor runs were carried out at the same j3 value on a 16^ x 48 lattice. This corre- 
sponds to the lattice spacing a = 0. 108(2) fm. We took two strange quark masses. For each of them, 
five up and down quark masses are taken covering the similar quark mass region. In total, ten runs 
of each length 2,500 HMC trajectories have been accumulated at 2 = 0. 



3. Topology issues 



If the topological charge is frozen in the QCD simulation, is there any problem? The answer to 
this question is obviously "yes". The QCD vacuum is required to be the 0-vacuum, a superposition 
of different topological sectors, in order to satisfy the cluster decomposition property. Therefore, 
we cannot sample the correct QCD vacuum unless the topological charge fluctuates sufficiently. 
This is a serious problem for everyone doing dynamical QCD simulation aiming at approaching 
the continuum limit, since the topological charge would not change near the continuum limit irre- 
spective of which fermion formulation one employs. 

We believe that the solution (one of the solutions, at least) to this problem is to reconstruct 



the 0-vacuum physics from the fixed topology simulations [gOp. This is based on the observation 
that the effect of fixed topological charge is a finite volume effect of 0{l/V), which can be derived 
from a simple Fourier transform between the 6 vacuum and the fixed Q "vacuum". For instance, 
the partition function at a given topological charge Q is obtained using a saddle point expansion as 
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Figure 1: Pseudo-scalar density correlator (3.3) summed over spatial separation and plotted as a function 
of t. The data for Nf — 2 QCD at am — 0.025. Solid line shows a fit with a constant giving the topological 
susceptibility. 



where Xt is the topological susceptibility and C4 characterizes a deviation from the Gaussian dis- 
tribution at a finite volume. Therefore, if one can calculate Xt and C4 (and higher order terms if 
necessary), the 0-vacuum can be reconstructed by simply adding the different topological sectors 
as Z{d) = T^Q^Q^ The method to calculate these key parameters is discussed shortly. 

The relation between the 6 and Q vacuums can be extended to any physical quantities. For 
example, for a (CP-even) two-point function, one can express it in a given topological charge Q 
in terms of its counterpart in the d{=0) vacuum G(0) and its second and fourth derivatives G*^^^(0) 
andGW(O) as pO, 



Ge = G(0)+GP)(0)^ 



1 



XfV 2xfV 



up to higher order corrections in l/V. Therefore, provided that the 0-dependence of the quantity of 
interest is known, the 0-vacuum physics can be reconstructed. The 0-dependence can be obtained 
for the quantities that can be analysed using chiral perturbation theory; if it is not applicable one 
needs the data at several Q to extract G^^' (0) for instance. 

By applying the formula ( p.2| ) for a topological charge density correlator, or equivalently a 
correlator of flavor-singlet pseudo-scalar densities mP{x), one obtains a relation 

jim (mPWmP(0))e = {xt + ' ■■) + O^e^^n'^), (3.3) 

where (• • •)e stands for an expectation value in the given Q. The topological susceptibility Xt can 
be extracted from the asymptotic value of the correlator in (^). One expects a negative constant 
when 2 = 0, which is intuitively understood as follows. When one finds a positive topological 
charge density at the origin, it is more likely that a negative value is observed at a far distant point 
X, if the net topology is fixed to 2 = 0. The value will be suppressed when the volume is increased. 
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Figure 2: Topological susceptibility as a function of sea quark mass. The data for Nf = 2 (circles) and 
Nf = 2 + 1 (triangles) are shown. The solid curves represent a fit with the form of chiral effective theory. 



An example obtained in two-flavor QCD is shown in Figure |I] p2[]. We observe a clear plateau 
for all six different ensembles with different sea quark masses. At am = 0.050, we obtain a consis- 
tent results from the configurations with different topological charges 2 = 0, —2, and —4. 

Sea quark mass dependence of Xt is expected to behave as = £/ ( 1 + 1 //n^/ + 1 / ) at the 
leading order of chiral perturbation theory [23]. In two-flavor QCD with degenerate sea quark mass 
m, it becomes Xt = mL/2. Figure^ shows the results for 2- and 2-1-1 -flavor QCD The data show 
a clear linear dependence on the sea quark mass. The fits to two- or 2-1-1 -flavor formula yield the 
value of chiral condensate Z. After a renormalization, we obtain £(2 GeV) = [242(5)(10) MeV]^ 
{Nf=2) and E(2 GeV) = [240(5)(2) MeV]^ {Nf=2+l). 



4. Physics applications 

4.1 Chiral condensate 

For the lattice calculation of the chiral condensate, the exact chiral symmetry plays a crucial 
role. Without the exact symmetry the scalar density operator on the lattice (i/ai/a)''^" may mix with 
the identity operator under radiative corrections. This leads to a power (cubic) divergence in the 
matching onto the corresponding continuum operator {yyY""' , which prevents one from calculat- 
ing this quantity on the lattice with any useful precision. 

Thanks to the exact chiral symmetry, we are able to precisely calculate this quantity from 
several different sources. They contain the use of the Banks-Casher relation (as demonstrated in 
Figure ^, the spectrum of low-lying eigenvalues (with a help of the chiral random matrix theory), 
the meson correlators in the £-regime, the topological susceptibility, and the Gell-Mann-Oaks- 
Renner (GMOR) relation. 

In the e-regime, where pion Compton wavelength is longer than the spatial extent of the lattice, 
the dynamical degrees of freedom is dominated by the zero-momentum modes of pions. In this 
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Figure 3: Spectral density p{X) as a function of X, calculated from lattice data in the e-regime. Banks- 
Casher relation infers a constant at A = 0, which is not of course satisfied at a finite volume, but the remnant 
is clearly seen. The horizontal dashed line is an expectation for a nominal value of E. 
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Figure 4: Chiral condensate calculated in two-flavor QCD with several different methods: e-regime eigen- 
value spectrum |5||, p-regime eigenvalue spectrum [p5||, e-regime meson correlator [p6[], topological 
susceptibiUty GMOR relation Q. 



circumstance, the low-lying eigeimiode spectrum can be obtained using the Chiral Random Matrix 
Theory (ChRMT). The chiral condensate can be extracted by matching the spectrum calculated on 
the lattice with the ChRMT prediction. Weobtainr^(2 GeV) = [251(7)(11) MeV]^ P3,|25l]. The 



meson correlator calculated in the £-regime can also be used to determine the chiral condensate 
I as well as the pion decay constant F. We obtain l'^^(2 GeV) = [240(4)(7) MeV]^ and F = 



87(6)(8) MeV [|26p. 

Figure ^ compares the chiral condensate extracted with various observables in two-flavor 
QCD. They are in good agreement. This strongly indicates that the low-energy dynamics of QCD 
is indeed determined by the pion degrees of freedom and therefore that the chiral symmetry is 
spontaneously broken, which is usually assumed when constructing the chiral effective theory. 
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4.2 and fn 

Near the massless limit of quarks the chiral effective theory is expected to provide a good 
description of the low-energy QCD. For larger pion masses, one must calculate loop corrections 
and also introduce irrelevant operators to subtract associated divergences. This leads to the chiral 
perturbation theory (ChPT), which is organized as an expansion in terms of small m| and p^. The 
region of the convergence of this chiral expansion is not known a priori. Using lattice QCD, one can 
test the expansion and identify the region of convergence. With the exact chiral symmetry, the test 
is conceptually cleaner, since no additional term to describe the violation of chiral symmetry has 
to be introduced. (With other fermion formulations, this is not the case. The unknown correction 
terms are often simply ignored.) 

For the pion mass rrijc and decay constant /;;: the expansion is given as 



2 

=2B 



1 2 

1 + ^x\nx + c-iX + 0{x ) 



(4.1) 

= f[l-x\nx + c^x + 0{x^)] , (4.2) 



where m-n and denote the quantities after the corrections while m and / are them at the leading 
order. The expansions ( p~l| ) and ( p3| ) may be written in terms of either x = 2m^/(4;r/)^, x = 
2ml/{4Kf)^, or = 2m\/{Anfnf- (we use a notation of fn = 131 MeV). Theoretically, they all 
give an equivalent description at this order, while the convergence behavior may depend on the 
expansion parameter. 

Figure ^ shows the comparison of different expansion parameters [27]. The fit curves are 



obtained by fitting the three lightest data points with the three different fit parameters. They all 
provide equally precise description of the data in the region of the fit. If we look at the heavier 
quark mass region, however, it is clear that only the -expansion gives a reasonable function and 
others miss the data points largely. This clearly demonstrates that at least for these quantities the 
convergence of the chiral expansion is much improved by the use of the -parameter than the 
conventional choice, i.e. the x-expansion. 
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Figure 6: Low energy constants of two-flavor QCD as determined from the (N)NLO analysis of our lattice 
data. The results are obtained using 4, 5, or 6 lightest sea quark mass data. Phenomenological numbers are 
taken from Colangelo-Durr p^, Colangelo et al. [^; the previous JLQCD results were obtained from the 
meson correlators in the e-regime [p6|]. 

Therefore, we consider only the ^ -expansion when we extend the analysis to the next-to-next- 
to-Ieading order (NNLO). I do not write down the NNLO terms here, but an important point is that 
the coefficient of the leading log terms of the form (i^ Ini^)^ is determined solely from ChPT and 
thus is not a fit parameter. The terms of the form and E, ^ have to be determined by the fit. 

It is indeed possible by a combined fit of m\/mq and Z^; we obtain the low energy constants as 
shown in Figure ^. 

It is clear from Figure ^ that the results depend on the order of the chiral expansion, i.e. NLO 
or NNLO, especially for the NLO low energy constants £3 and £4. This is because the NNLO term 
appears with a relatively large numerical coefficient, such as ^31^(1 — 9/2t, Int,), for instance. A 
large effect is then expected for the determination of unless is much less than 0.1. The NNLO 
terms are mandatory in the analysis of lattice data. 

We are currently extending the NNLO analysis to the partially quenched data in two-flavor 
QCD and to the 2-1-1 -flavor QCD data [§. 

4.3 Pion form factor 

Pion vector and scalar form factors provide another simple testing ground for the chiral dy- 
namics of pion. They are defined as 

{K{p')\V^\K{p)) = {p^+p'^)Fv{q^), (4.3) 
{7t{p')\S\7tip)) = Fs{q^), q^^p'^ - p^, (4.4) 

for the vector current and the scalar density S. At small momentum transfer q^, the form factors 
are characterized by the charge and scalar radii, {r^)y, and (r^)^, as 

M^') = l + Ur')W + Oiq''), (4.5) 
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Figure 7: Pion vector (left) and scalar (right)form factors calculated at am ~ 0.050. 
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Figure 8: Vector (left) and scalar (right) radii of pion as a function of M^. Solid curves show the fits with 
the NNLO ChPT formulae. Their breakdown into NLO (dashed) and NNLO (short dashed) contributions 
are also shown. 



Fsiq^) = Fs{0) 



(4.6) 



Note that the vector form factor Fv{q^) is normalized at = due to the current conservation. 

We use the all-to-all quark propagator technique to calculate the form factors ^. The 
main advantage of this technique is, of course, that one can calculate the disconnected diagram 
contribution to the scalar form factor. Furthermore, it substantially reduces the statistical error in 
the calculation through an average over the source point. 

We calculate the form factors at four lightest sea quark mass configurations available in two- 
flavor QCD. The data at a fixed quark mass are shown in Figure ^ For the vector form factor we 

attempt a fit ansatz motivated by the pole dominance Fv{q^) = 1/(1 — /niy) +ciq^ H with 

niy the vector meson mass obtained at the same quark mass. This is justified by the analyticity 
of the form factor and also by a reasonable assumption that the higher pole contributions are well 
approximated by a Taylor series. The data are, in fact, nicely fitted with this function as shown in 
Figure ^ (left). The higher pole effects are small but visible. For the scalar form factor there is no 
obvious pole, and we use a simple polynomial of q^ to fit the data (right panel). 

The charge and scalar radii thus extracted are plotted as a function of m| in Figure ^ together 
with the fits with the NNLO ChPT formulae. For these quantities, the ChPT predicts a logarithmic 
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divergence in the chiral limit ^ (l//|)lnm^ with a fixed numerical coefficient. The data do not 
show the expected divergent behavior, but the fit is possible and gives a value consistent with 
the phenomenological value at the point of physical pion mass. This suggests that the divergent 
behavior shows up in the even smaller quark mass region. We note that the fit is possible only when 
we include the NNLO terms, which contain additional logarithmic divergences in the chiral limit. 

4.4 Vacuum polarization functions 

The vacuum polarization functions defined through 



(4.7) 



contain rich information of QCD dynamics in both perturbative and non-perturbative regimes. (J^ 
is either a vector or an axial-vector current. We restrict ourselves in the flavor non-singlet 
currents, though the flavor indices are suppressed in the expressions.) Experimentally, the vacuum 
polarization function can be extracted from e^e^ annihilation or from T decay processes. 

The spontaneous chiral symmetry breaking can be probed with the vacuum polarization func- 



tions using the Weinberg sum rules [33], such as 



Jit 



lim Q 
lim - — 



n 



(1+0) 



n 



1+0) 



n 



(1+0) 



(4.8) 
(4.9) 



where S stands for the Peskin-Takeuchi's parameter used in the context of the precision electroweak 



measurement [)34|]. In chiral effective theory, it is related to a low energy constant Liq. These sum 
rules can be obtained in the limit of massless quarks. Because of their VV —AA structure, they 
vanish when the chiral symmetry is not spontaneously broken. 

Along these lines, another interesting quantity is the electromagnetic mass difference of pion, 
which is expressed in this Umit as 



Ami 



3 a, 



'EM 



n 



(1+0) 



(4.10) 



which is known as the Das-Guralnik-Mathur-Low- Young sum rule [35]. 

For a lattice calculation of these quantities, exact chiral symmetry is essential, since we have 
to extract a tiny difference between the vector and axial-vector channels. In the calculation using 
the overlap fermions, we checked that even the lattice artifacts precisely cancel between V and A 



In Figure ^ we plot the results for the difference 11^^ {Q 



The chiral limit is taken 



by assuming a linear quark mass dependence (solid curve) or by using the one-loop ChPT formula 
(dashed curve). The use of ChPT is slightly questionable, because the lowest (non-zero) value of 
on our lattice is already as large as (320 MeV)^ and the second lowest ~(650 MeV)^ is clearly 
out of the range. The high region, on the other hand, may be analyzed using the operator product 
expansion (OPE). OPE is also used to extend the region of to infinity as required in ( [4.10| ). For 
the pion mass difference, we finally obtain m| = 992(12)(+j3° )(149) MeV^ (see [g§] for details), 
which may be compared with the experimental value 1,261 MeV^. 
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Figure 9: Difference of the vacuum polarization functions TIy\q^) — Ti''J^\q^) muhipHed by Q^. Data 
are shown for different quark masses; the chiral hmit is also plotted by a solid (linear extrapolation) and a 
dashed (chiral logs) curve. 
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Figure 10: Vacuum polarization function Tly^^ {Q^) as a function of Q^. The lattice data at different quark 
masses are shown together with the fit curve with the OPE formula (solid curve). 



Another use of the vacuum polarization function is an extraction of the strong coupling con- 
stant by matching the lattice data with the OPE expression in the perturbative regime. For this 
purpose we consider the Adler function Dj{Q^) = —Q^dYlj{Q^)/dQ^, which is free from ultra- 
violet divergence, and thus renormalization scheme independent. It means that one can use the 
expressions obtained in the continuum perturbation theory to describe the lattice data. (In the prac- 
tical analysis, we directly use Ylj{Q^) instead, to avoid numerical derivative in terms of Q^. We 
then float a divergent constant as a fit parameter.) 

In the lattice calculation, it is necessary to subtract lattice artifacts mainly due to non-conserved 
lattice currents. Non-perturbative subtraction is possible by utilizing several different momentum 
configurations giving the same [|^. 

Lattice results are shown in Figure IC. The data in high regime are fitted with the pertur- 
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bative expression supplemented by OPE 



where perturbative expansion of Co, Cm, etc is known to a} in the MS scheme, c is the divergent 
constant term. 

In the analysis, we may use the value of the chiral condensate {qq) obtained through other 
quantities as described above. We then obtain = 234(9)(+'o) MeV in two-flavor QCD 
This result is compatible with previous lattice calculations in two-flavor QCD, 250(16)(16) MeV 



[38] or 249(16X25) MeV [39]. 



4.5 Nucleoli structure 

The last application of our dynamical overlap fermion simulations reported in this talk is a 
calculation of the nucleon sigma term 

OnN = mud{N\uu + dd\N) (4.12) 

and the strange quark content 

, = ^m^. (4 13) 

{N\uu + dd\N) 

For GjcN both connected and disconnected quark lines contribute, while only the disconnected dia- 
gram contributes to the strange quark content. 



Instead of directly calculating the matrix elements in (1.12) and (4.13), we utilize the Feynman- 
Hellman theorem. Namely, we calculate the derivatives of nucleon mass in terms of valence or sea 
quark mass. They are related to the connected and disconnected contributions as 



5 '"sea 



{N\uU + dd\N) conn, (4.14) 
{N\au + dd\N)di,c (~ 2{N\ss\N) ) , (4. 15) 



where the subscripts "conn" and "disc" stand for the connected and disconnected contributions to 
the matrix element, respectively. The disconnected piece becomes identical to the strange matrix 
element 2{N\ss\N) when the sea quark mass is set equal to the strange quark mass. To obtain 
the derivatives we need partially quenched (myai 7^ '«sea) data set. For the fit of the quark mass 
dependence, we use the partially quenched chiral perturbation theory calculated at one-loop for the 

nucleon mass ©Ell]. 



Figure |n| (left panel) shows the nucleon mass calculated on our two-flavor QCD configura- 
tions as a function of m| [42]. Finite volume effect is corrected using the expectation from ChPT 
The data are nicely fitted with the ChPT formula at 0{p^) or at 0{p^), through which we 



obtain o„n = 52(2)(;20)(^+5) j^^y ^^^^ further details, see \^). 

The derivatives can be obtained after fitting the partially quenched data set with the corre- 
sponding ChPT formulae. The results are shown in Figure [n| (right panel). The plot clearly shows 
that the disconnected contribution is smaller than the connected contribution. In particular, the dis- 
connected piece is tiny for the mass corresponding to the strange quark. As a result the y parameter 
we obtained is smaU: y = 0.030(16)(;^)(+^). 
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Figure 11: (Left) Nucleon mass in two-flavor QCD. The data after the finite volume correction are shown 
by red symbols. The curves represent the fit with baryon ChPT formulae. (Right) The connected (upper) 
and disconnected (lower) contributions to the sigma term as a function of the quark mass. Sea and valence 
quarks are taken to be equal. Physical up/down and strange quark masses are shown by vertical lines. 



Previous lattice calculations gave rather large values fory: 0.66(15) [44], 0.36(3) [45], 0.59(13) 
[%]. The problem in these calculations is that the critical mass depend significantly on the sea 
quark mass when one uses the Wilson-type fermion due to the lack of chiral symmetry. One can 
misidentify this unphysical effect with the physical sigma term. (We note that the problem persists 
even when one directly calculates the matrix element rather than using the Feynman-Hellman the- 
orem. That is because the scalar density operator made up with valence quarks may mix with that 
with sea quarks. One has to subtract this operator mixing to obtain the correct results [|5|].) 

After subtracting this effect, UKQCD obtained a value consistent with zero, y = —0.28(33), 
albeit a large error due to the subtraction [4-7]. Our setup is free from this problem because of the 
use of the overlap fermion for both valence and sea quarks. The exact chiral symmetry plays a key 
role here, too. 



5. Conclusions 

In this talk 1 demonstrate that the dynamical overlap fermion may provide a clean approach 
to the problems related to the chiral symmetry and its spontaneous breaking. Despite the large 
numerical cost it requires, large scale simulation is feasible with the computational resources of 
0(10 TFlops). So far, we have completed runs on a 16^x48 and started test runs on a larger 
volume, 24^x48. 

The key to this success is the new strategy of treating the global topological charge, i.e. we fix 
the topology during the simulation and reproduce the 0-vacuum physics by correcting an induced 
finite volume effect. In doing so, the key observation is that the topological susceptibility can be 
correctly extracted from simulations at a fixed topology. 

The problem of topology freezing will become common among all lattice simulations as one 
sufficiently approaches the continuum limit, simply because this is the property of the continuum 
QCD. Our work is a first attempt to overcome that situation. 
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The exact chiral symmetry plays a crucial role in the calculation of many important physical 
quantities. For instance, we have calculated the chiral condensate using several methods for the 
first time and confirmed the consistency among them, which provides the most fundamental test of 
the chiral effective theory. We are extending the test to various quantities, such as pion mass, decay 
constant, charge radius and so on. (A calculation of the kaon B parameter is already published 



[ftSQ.) There are other interesting applications for which the exact chiral symmetry is required; so 
far we studied the VV — AA vacuum polarization and the nucleon sigma term. We plan to further 
pursue in this direction. 

The author is supported in part by Grant-in-aid for Scientific Research (No. 18340075). 
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